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Abstract 

If two parties, Alice and Bob, share some number, n, of partially 
entangled pairs of qubits, then it is possible for them to concentrate 
these pairs into some smaller number of maximally entangled states. 
We present a simplified version of the algorithm for such entanglement 
concentration, and we describe efficient networks for implementing 
these operations. 

1 Introduction 

The state of a single pure quantum bit, or qubit, is described by a vector in a 
2-dimensional Hilbert space spanned by basis vectors |0) and The state 
of n pure qubits (i.e. an n-qubit register) is described by a vector in a 2"- 
dimensional Hilbert space which is the tensor product of the 2-dimensional 
spaces for the states of each of the n qubits. Consider a 2-qubit register 
in a state described by the vector = -^lOO) + We call a pair 

of particles in this state an EPR pair, named after Einstein, Podolsky and 
Rosen, who discussed such particle pairs in their 1935 paper |[EPR35|| . It can 
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easily be shown that this vector cannot be factored into a tensor product of 
two 1-qubit states. That is 

|vl>) = -L|00) + i=|ll) ^ (ao|0) + ai|l)) ® (6o|0) + 

for any ao, ai,bo,bi. The amount of entanglement present in a bipartite 
quantum system can be quantified, and for this purpose we will treat a single 
EPR pair as possessing one unit of entanglement. 

In many scenarios involving quantum communication, an essential ingre- 
dient is the sharing of an EPR pair by Alice (the 'sender' of some information) 
and Bob (the 'receiver'). For example, when Alice and Bob share an EPR 
pair, they are able to perform quantum teleportation, a process useful for 
communicating quantum information. Using protocols involving the shar- 
ing of EPR pairs, some distributed computation tasks can be achieved using 
fewer bits than could be achieved using only a classical channel (see e.g. 
| BCW98|| and |P^az99|| ). Suppose Alice and Bob share a known entangled 



pair of qubits 

1^) = aoo|0)|0) + aoi|0)|l) + aio|l)|0) + an|l)|l), 

where the first qubit is in Alice's possession and the second qubit in Bob's. 
The Schmidt decomposition for this bipartite system allows us to express the 
state of this pair of qubits as 

1^) = a\ao)\bo) + p\ai)\bi), 

for some non-zero positive real numbers a and /3, and unit vectors |ao) and 
\ai) that form a basis for Alice's system, and unit vectors |foo) and that 
form a basis for Bob's system. Since Alice and Bob can each locally perform 
the one-qubit unitary operations 

|«o) ^ |0),|ai) ^ |1) 

and 

|f>o)^|0),|6i)^|l) 

respectively, we will assume that Alice and Bob share an entangled state of 
the form 

a|00) 
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If |a| = 1/3 1 = then the state is an EPR state, and is said to be 
maximally entangled. If \a\ 7^ \j3\ then the state is less entangled, and if 
either |a| or \P\ equal 0, then the state is completely non-entangled. 

Consider a 2n-qubit system of the form |\&) = (a|00) + shared 
by two parties, Alice and Bob, where |a| 7^ |/3|. Now suppose Alice and 
Bob want to share some maximally entangled EPR pairs for some commu- 
nication task. A natural question is "how many EPR pairs can Alice and 
Bob distill out of |\E'), performing local operations and communicating clas- 
sically" ? An upper bound on the expected number of EPR pairs that can be 
distilled is the "entropy of entanglement" of |\E') defined to be von Neumann 
entropy of either pA = Ttb\'^){'^\ or pB = Tr^|\I')(^|. These quantities are 
both equal to the Shannon entropy of the eigenvalues of (|\E') (\E'|)'^ (which 
are the squares of Schmidt coefficients of the state I^E')'*). This quantity 
equals n times the von Neumann entropy of |\E')(\E'|, namely nH{\a'^\), where 

H{p) = plog + (1 — p) log (^z^^ ■ Foi' example, the Von Neumann en- 
tropy of an EPR pair is H{\{^f\) = 1. 

The process of distilling EPR pairs out of |\E') is called entanglement 
concentration. Local operations for performing entanglement concentration 



have been by Bennett, Bernstein, Popescu and Schumacher in IPBPS95 
The expected amount of concentrated entropy of entanglement is 



x:i«^i"-(i-i«^iy(")iog, 



(1) 



and they show that this quantity is in nH{\a'^\) — 0(log?7,). 

In section § we describe the approach detailed in |PBPS95[] . We then 
describe a new way of extracting a specific number of EPR pairs instead of 
the method suggested in |[BBPS95"| . In section ^ we will give a description 
of a quantum network for performing the main local basis change necessary 
for performing entanglement concentration. In section ^ we summarise how 
to implement entanglement concentration. 
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2 Local Operations for Entanglement Con- 
centration 

Consider concentrating the entanglement of the state |^) defined above, and 
without loss of generality we assume that q;,/3 are positive real numbers. 
Consider the case for n = 3 qubits: 



Separating Alice's qubits from Bob's, we can re-write the above state as: 



a^lOOO) |000) + a=^/3(|001)|001) + |010)|010) + |100)|100)) 
+a/?2 (|011)|011) + |101)|101) + |100)|100)) + 

In general, if we have n copies of q;|00)-|-/3|11), by appropriately reordering 
the qubits we get: 



where Ahce's qubits are labelled with an "a" and Bob's are labelled with a 
"6" , and H(x) is the number of Is in the string x, also known as the Hamming 
weight of X. On the right hand side of the equality, the state is written in 



{a\oo) + p\n)f 



= 0)1 0)1 0) 
W/?(|00)|00)|11) + |00)|11)|00) + |11)|00)|00)) 
+aP^ (|00)|11)|11) + |11)|00)|11) + |11)|11)|00)) 
+/3=^|11)|11)|11). 



alice bob alice bob alice bob 



alice 



bob 
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terms of the symmetric basis. The symmetric space is an (n+ l)-dimensional 
subspace of the 2"^-dimensional state-space for the register. The i^^ symmetric 
basis state is a uniform superposition of the computational basis states having 
Hamming weight i. Ahce and Bob can each measure the Hamming weight of 
their half of the state The measurement is implemented by introducing 

an ancilla of size C'(logn). A sequence of controlled- [add 1] operations is used 
to add the Hamming weight of each qubit of |\E') into the ancilla. This is 
implemented by network shown in Figure |I|. 




|o) — 
|o> — 



Figure 1: Network to compute the Hamming weight. \xi - ■ ■ x„) |00 ■ ■ ■ 0) 
\xi ■ ■ ■x„)|H(a;i ■ ■ -x^)). 



+ 1 



+1 



+ 1 



|//(Xi..j:„)) 



Suppose Alice measures the Hamming weight of |\E') and obtains the result 
\j) (Bob will measure the same j whenever he performs the same measure- 
ment). This state after the measurement is 




which can be thought of as a superposition of (") n-bit strings. (Of course, 
the measurement is not necessary, and the remainder of the algorithm could 
be controlled quantumly upon the value j.) Let r = [log2 (")]. Define a 
function / on these (") strings that maps the strings of length n with 



5 



Hamming weight j (in lexicographic order) to the integers from to — 1: 

/(OO . . . 00n_^) = 00 ... 

j 

/(OO . . . ion_^) = 00 ... 1 



/(11_^00...0) 

j 



Wc can extend / so that it defines a permutation of all n-bit strings. Then 
we have: 
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ix)|x) M j2 



ij) H(a;)=j H(a 

= E io) ly)lo)ly)- 

" n—r r 

If (p = T' , then ignoring the first n — k bits on both sides gives us 

2''-l 

y=o 

which is r EPR-pairs, and the entanglement of |^) has been concentrated. 

However, in general, (") will not be a power of 2. Let A;= [logs (")J + 1. 
We describe a quantum network that will produce some number Q <l <k — l 
of EPR-pairs (we use this definition for k in place of the previous definition for 
r for convenience in describing a network that will behave the same whether 
or not (") is a power of 2). The expected number of EPR pairs will be at 
least k — 2. We illustrate this for n = 3 entangled pairs of qubits. Consider 
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the binary representation (") = X2X1X0 = X2 ■ 2^ + xi ■ 2^ + xq ■ 2°. We have 

y=0 y=0 y=X2-22 y=X2-22+xi-2i 

(2) 

Notice that if a;2 = 1 then the above sum includes 000 < y < Oil. These 
are included in ^y=o ~^ |y)|y) which is the first term on the right side of (H) 
(if X2 = 0, then this term is empty). Similarly, if xi = 1 the sum includes 
X2OO < y < X2OI and if xo = 1 it includes X2X1O < y < X2X1I. So we can 
write the sum (Q) as follows: 



1 



oil a;201 X2X1O 



E ly)ly) = ^2 E + E ly)ly) + ^o E 1^)1^)- (3) 

y=0 y=000 y=X200 y=a;2a;i0 

In other words, for each j such that Xj = 1, we have the superposition of 2^ 
strings. Alice and Bob wish to project to one of these superpositions of 2^ 
strings, since that will provide them with j EPR pairs. 

The first term on the right side of @ contains the strings 000, 001, 010, 
Oil; all the strings beginning with a 0. Suppose Alice performs a measure- 
ment of the qubit in the leftmost position (i.e. corresponding to 2/2) of her 
(")-i 

share of the state ^y=o |y)|y) (Bob will obtain the same result whenever 
he performs the analogous measurement on his share). In addition, Alice also 
has the corresponding bit X2 in a register containing the binary expansion of 
(") . There are three cases to consider: 

CASE 1: y2 = and X2 = 1. 

In this case the joint state after the measurement is 

oil 

E 

y=000 

which is the first term on the right side of (|^). Ignoring 1^2), this is 2 
EPR pairs. 

CASE 2: 2/2 = and X2 = 0. 
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The state after the measurement is 



001 OxiO 
y=000 y=Oa;iO 

Ignoring the leftmost qubit \y2) this is equal to 

01 xiO 
y=00 y=a:iO 

CASE 3: y2 = 1. 

In this case we know X2 = 1. So the post- measurement state is 

101 la:iO 

xi J2 ly)ly) + ^o J2 ly)ly)- 

y=100 y=la;iO 

Ignoring the leftmost qubit \y2) this is equal to 

01 xiO 

xi Y |y)|y) + ^0 Y 

y=00 y=a:i0 

If Alice's measurement results in case 1, then the entanglement has been 
concentrated, and she makes no further measurements. Cases 2 and 3 both 
leave Alice and Bob with the state xi EJLoo ly)|y) + Zly=xiO |y)|y)- 1^ 
either of these cases, Alice discards the leftmost qubit \y2). She then repeats 
the measurement procedure, where this time the leftmost bits being measured 
are yi and Xi. The analogous three cases are considered again. 

This time case 1 would result in the post measurement state Y^y=oQ |y) |y)- 
Ignoring the leftmost qubit yi, this gives 1 EPR pair, and the procedure stops. 
Cases 2 and 3 both result in the post measurement state X^y^Q |y) |y), 
giving EPR pairs. 

It is easy to generalise this approach for k = [log2 (")J + 1. Alice (or 
Bob) measures (locally) the qubits yk-i, ■ ■ ■ ,yi, from "left-to-right", at each 
step checking the value of the corresponding bit Xi in the binary expansion 
of ("). She does this until, at some iteration I (where the first iteration is 
indexed 0), she finds = |0) and the corresponding bit x^-i-i = 1. 

When this occurs the procedure stops, having distilled k — I — 1 EPR pairs. 
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Figure 2: Network to measure how many EPR pairs have been distilled. 



A quantum network implementing the procedure is shown in Figure 
Since the Hamming weight j has been measured, (") can be efficiently com- 
puted. The binary representation of (") is encoded in a register \xk-i ■ ■ ■ Xq). 
The network makes use of an ancilla of size k, initially in the state We 
refer to this ancilla as the "control ancilla", and label its qubits of the by \ti) 
for < i < k — 1. For each i, \ti) is switched to |0) if both \yi) = |0) and 
\xi) = This is achieved using a sequence of doubly controlled NOT gates 
in the first stage of the network, where the NOT is applied to the target 
qubit if the first control qubit is in state |1) and the second control qubit 
is in state |0). Another ancilla of size 0{\og{k — 1)), which we will call the 
"measurement ancilla" is initially in the state \k — 1). In the second stage 
of the network, the value of the measurement ancilla is decremented by a 
sequence of controlled- [subtract 1] gates, controlled successively on each of 
the \ti) in the control ancilla. The net effect of the first two stages of the net- 
work is to decrement the measurement ancilla by one for each pair \yi)) 
until one such pair is found with = \l),\yi) = |0)). After such a pair is 
encountered, the measurement ancilla is not decremented any more. In order 
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to reverse the effect of any coupling that the network may have introduced 
between the primary register |y) and the ancilla |t), the same sequence of 
doubly controlled NOT gates that was used in the first stage of the network 
is applied again in the third stage. 

The control ancilla has been reset to its initial state by the third stage 
of the network, and the register |x) containing the binary expansion of 
is in a fixed computational basis state, since the value of j was fixed by 
the Hamming weight measurement performed earlier. Ignoring the state the 
control ancilla and the register |x), the joint state of Alice's system |y) and 
the measurement ancilla just before the final measurement is 

Xfc.i 5^ |y)|A;-l) + Xfc_2 Yl \y)\k-2) + --- + xo Yl 1^)10)- 

y=0'' y=Xk-iO''-'^ y=Xk-i---xiO 

The string Xk-iXk-2 ■ ■ ■ Xk-i correspond to the the leftmost / bits in the 
binary representation of (") . After the measurement of the control ancilla in 
the computational basis, the state is 

a;fc_ia;fc_2...a:fc_i0l'="'-l 

E ly)l^-^-i) 

y=Xk--i_Xk-2-Xk-fi''~'' 

for some < / < A; — 1. Ignoring the leftmost / + 1 qubits, this is 

-^k-l-l 

E ly)l^-^-i)- 

y={)k-l-l 

Each ignoring their respective leftmost / + 1 qubits, the joint Alice-Bob state 
is 

]k-l-l 

E \y)\y) 

y=Qk-l-l 

which is A; — / — 1 EPR pairs. Note that the state of Alice's ancilla (and 
Bob's, if he performs the same measurement procedure on his share of the 
state) indicates the number of EPR pairs that have been distilled. 

It should be noted that Alice and Bob can each carry out the above 
procedure locally, and they will obtain the same results. Alternatively, Alice 
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could perform the Hamming weight computation locally and send the result 
to Bob. Alice and Bob would both perform the permutation / by the method 
detailed in section In the last stage of the procedure, either Alice or Bob 
could perform the computation to determine the number of EPR pairs that 
have been distilled, and send the result to the other. 
It can be shown that given the superposition 

y=o 

the average number of EPR-pairs produced using this approach is: 

3 Implementing the permutation / 

The key step in the entanglement concentration protocol is the permutation 
/. We need to know how to implement this function. Recall that we start 
with a superposition of (") strings x, each having j Is and n — j Os. We 
want / to impose a lexicographic ordering on these strings. 
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Consider the following: 

00 ... 



00 ... 00 1 1 ... 1 



00 ... 1 



... 10.11 ... 1 

J-1 



n — 1 



n-\\ f- 
3 



,1 1 ... 1 ... 

3 




- 1 



/-I 




0. 



The first . ^) strings have a in the first bit position, and the remaining 

strings have a 1 in the first bit position. Define [y]nj to be the y*'^ largest 
string (treating the string as an integer represented in binary) of length n 
with Hamming weight j. Using this notation, the method for implementing 
/ is captured by the following recurrence: 



[y]n,i = 0[y]n-l,i 



if < y < 



n 



= 1 



n — 1 
3 



n-l,j-l 



n 



.fl"-^<y< 

3 J \3 



(4) 



We describe how the permutation can be implemented on a quantum 
computer. Let [y]n,j = ^162 . . . Then start with a string between 00 ... 
and (") — 1, and ancilla holding the values n, j, and a space for the output 
strings /"^(y): 



|y)|n)|j)|00...0). 
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Apply an operator T which performs the following mapping: 

T I 71 — 1 

|y)|n)|j)|00...0)^|y)|n)|j)|00...0) if < y < ( ^ 



T 



|y)|n)|j)|10...0) if ( ^ 



n — 1\ (n 



The result is: 



|y)|n)|j)|00...0) ^ |y)|n)|j)|M...O) 



for some h\ G {0,1}. Then perform the following subtraction operation S', 
controlled quantumly on the value of h\. 

\y)\n)\3) ^ \y)\n - \)\3) if 6i = 

-^|y-(''T^))|n-i)b-i) if 6i = l. 

Then repeat T and iS, this time on only the right-most n — 1 bits of the 
registers \y) and |6iO..O). Applying T and S in this way, a total of n times, 
realises the recurrence (1), and gives us an implementation of 

|y) ^ l[y]nj) 

and thus the same network maps 

(•;)-■ 

E iy> - E 

y=0 /^(x)=i 

The permutation / is realised simply by running this procedure backwards. 



4 An algorithm for entanglement concentra- 
tion 

We now have the tools to state an algorithm for implementing the entangle- 
ment concentration protocol described in section 4.1. The algorithm is the 
following: 
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1. Begin with the state |^) = (a|00) + 

2. Ahce and Bob each perform a Hamming- weight measurement on their 
half of obtaining the same result 

3. Alice and Bob each perform the permutation / on the resulting super- 
position. 

4. Alice and Bob each use the network of Figure ^ to determine how many 
EPR pairs they share. 

5. The result is some known number of perfect EPR pairs. For a par- 
ticular j, the expected number is between k — 2 and k — 1 where 
k=[\og, (pJ+1. 

Each of the above steps have been detailed in the preceding sections, and so 
we have a complete description of the implementation. Since the probability 
of measuring \j) in step || is 

i«r-(i-i«^ir(")> 

the expected number of EPR pairs is at least 

|:i«T-'(i-i«v(;)(i.o^;)j-i) 

and comparing to equation (||) shows that we get at least 

nH{\a^\) - 0{\ogn) 
EPR pairs on average. Note that the theoretical maximum is nH{\a'^\). 
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